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ABSTRACT 

Observational biases against finding Milky Way (MW) dwarf galaxies at low Galactic latitudes (b 
£ 20°) and at low surface brightnesses (/iy,o i> 29 mag arcsec -2 ) currently limit our understanding of 
the faintest limits of the galaxy luminosity function. This paper is a proof-of-concept that groups of 
two or more RR Lyrae stars could reveal MW dwarf galaxies at d > 50 kpc in these unmined regions of 
parameter space, with only modest contamination from interloper groups when large halo structures 
are excluded. For example, a friends-of-friends (FOF) search with a 2D linking length of 500 pc could 
reveal dwarf galaxies more luminous than My = -3.2 mag and with surface brightnesses as faint as 
31 mag arcsec -2 (or even fainter, depending on RR Lyrae specific frequency). Although existing 
public RR Lyrae catalogs are highly incomplete at d > 50 kpc and/or include <1% of the MW halo’s 
volume, a FOF search reveals two known dwarfs (Bootes I and Sextans) and two dwarf candidate 
groups possibly worthy of follow-up. PanSTARRS 1 (PS1) may catalog RR Lyrae to 100 kpc (in the 
absence of Galactic extinction) which would include up to ~15% of predicted MW dwarf galaxies. 
Groups of PS1 RR Lyrae should therefore reveal very low surface brightness and low Galactic latitude 
dwarfs within its footprint, if they exist. With sensitivity to RR Lyrae to d /> 600 kpc, LSST is the 
only planned survey that will be both wide-field and deep enough to use RR Lyrae to definitively 
measure the Milky Way’s dwarf galaxy census to extremely low surface brightnesses, and through the 
Galactic plane. 

Subject headings: galaxies: star clusters — galaxies: dwarf — stars: distances — stars: variables: 
other — techniques: photometric 


1. INTRODUCTION 

The Milky Way’s (MW’s) dwarf galaxy population pro¬ 
vides a unique laboratory to study galaxies at the bot¬ 
tom of the galaxy luminosity function, and to quan¬ 
tify tensions between cosmological predictions and ob¬ 
servations on sub-galactic scales. For example, the 
discrepancy between the thousands of dwarf galaxy- 
mass dark matter halos p redi c ted to orbit the Milky 
Way ( e.g. iSpringel et al.1120081 iGarrison-Kimmel et al.1 
l2014al) and the ^35 dark matter do minated dwarf 


galaxies observed to orbit the MW (IMcConnachie 

2012 

: The DES Collaboration et al.ll2015HKonosov et al. 

2015 

) has been long referred to as the “missing satcl- 


lite problem” dKlvnin et all 119991 iMoore et al.Ml999ll . 


The discrepancy between the predicted and observed cir¬ 
cular velocities of the most massive MW dwarf satel¬ 
lite^ has been coined the “too big to fail” problem 
(jBovlan-Kolchin et all 120111 120121 ) . Most recently, in¬ 
creasing attention has been paid to possible differences 
between the predicted and observed spatial distributions 
of dwarf galaxies ar ound the MW, M31 and other MW- 
analog galaxies (e . g.lKroupa et al . 2005; Pawlowski et al.1 
l2012t llbata et al.1l20l4 iPhillips et alJl2015l) . Resolving 
any of these open puzzles about the MW’s dwarf satel¬ 
lites, or measuring the empirical relationships needed to 
understand galaxy formation at its lower limit (for exam¬ 
ple, the empirical size-luminosity-distance relationship), 
requires a more complete and unbiased census of MW 
dwarfs. 

The current detection biases against low surface bright- 
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ness and low Galactic latitude dwarfs create an incom¬ 
plete picture of the MW’s dwarf galaxy population. For 
example, some theories predict th e existence of extre mely 
low surface brigh tness “stealth” (iBullock et al 1120 Kill and 
“fossil” galaxies dBovill fe Ricottill2011h thatwould have 
evaded detection with existing techniques. Figure |T| 
shows that no dwarfs are known to exist in the larger 
size - lower luminosity space that such extreme dwarf 
galaxies may occupy. This demographic gap may reflect 
something physical about the formation of low luminosity 
dwarf galaxies, or may simply reflect the detection bias 
against dwarfs with V-band centra l surface brig ht nesses 
fainter than ~29 mag arcsec -2 (e.glKoposov et al .1120081 : 
IWalsh et al.1120091 iThe DES Collaboration et al.112015T) . 
Moreover, MW dwarf galaxies are detected and stud¬ 
ied by their resolved stellar populations. The increasing 
number of foreground stars at Galactic latitudes lower 
than b ~ 20° therefore limits the detectability of dwarf 
galaxies over a si gnificant fraction of the MW’s volume 
(IWalsh et alJl2009l) . 

RR Lyrae (RRL) stars may provide a new means to 
detect previously unseen MW dwarf galaxies. RR Lyrae 
are short-period variable stars found in old and metal- 
poor stellar populations. At least one RR Lyrae has been 
found in each Milky Way dwarf galaxy searched for such 
stars (see Section 13.11) . There are two main categories 
of RR Lyrae stars: type c RR Lyrae stars (RRc), which 
pulsate in the first overtone and have periods between 
0.25 and 0.4 days, and type ab RR Lyrae stars (RRab), 
which pulsate in the fundame ntal mode and have peri¬ 
ods between 0.4 and 0 .8 days (|Smithl [T995fc iVivas et al.1 
12004 ISesar et, al .1120071) . RR Lyrae have light curves that 
exhibit distinct shapes, periods and amplitudes, making 
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Fig. 1.— The observed size-luminosity-distance relations of confirmed (closed points) and recently published candidate (open points) 
Milky Way dwarf galaxies. A line of constant central V-band surface brightness is overplotted on the Left panel. Observational bias against 
finding very low surface brightness dwarfs (upper left portion of the left-hand panel) and distant, faint dwarfs (bottom right portion of the 
right-hand panel) continue to limit our census of Milky Way dwarfs. This figure includes 33 MW satellites: Hercules, Boo I, Boo II, Leo 
IV, Leo V, Pisces II, CVn I, CVn II, UMa I, UMa II, Coma Berenices, Wil l 1, Seg 1, Seg 2, D raco, Fornax, Sculptor, Carina, Sextans, 


Le o I, Leo II, Ursa Minor JSand et a hm a nd references therein) , Leo T Jde^on^et^l 
2 dThe D ES C ollaboration et a lj|2015|) , Gru 1 JKoposov et al.ll201 “‘' 

JLaevens et *30201511 . Hydra 2 IIMartin et al.H2015l ). 


in), s===u= 

15), Hor 2 ( Kim V Jeri enl 


2005), Ret II, Eri 2, Tuc 2, Hor I, Pic 1, Phe 
2015 i. Peg 3 lIKim et aI.II2015f ). Triangulum 3 


them relatively easy to identify in time domain surveys. 
RRab stars are also particularly good standard candles, 
with a well calib rated l umi nosity -[Fe/H] relation at opti¬ 
cal magnitudes (lChaboveiill999l) : 

M v = (0.23 ± 0.04) [Fe/H] + (0.93 ±0.12) (1) 

These combined properties make RR Lyrae stars valu¬ 
able tracers of structure throughout the Milky Way’s 
halo. Until now, RR Lyrae have primarily been used 
to map streams and unbound overdensities in the Milky 


Way’s halo (e.g.lVivas & Zind 2006); 

Watkins et al. 

2009; 

Brake et al.1 2013al ISesar et aLl|201£ 

: Buffau et all 

2014; 

Simion et al. I2014J). However, as wide-held time domain 


surveys become sensitive to RR L yrae at increasin g dis¬ 
tances into the MW’s halo (e.g. lOluscvi ct al]|2012i ). RR 
Lyrae could be used as disc overy tool s fo r dwa rf galaxies 
themselves (see Section [2]). iSesar et ahl (12014D proposed 
that single RR Lyrae in the outer halo could be bea¬ 
cons of underlying, extremely low luminosity (L < 1000 
Lq) dwarf galaxies that are otherwise invisible. They 
demonstrate that stacked star catalogs around single, 
distant RR Lyrae stars provide interesting constraints on 
the presence of extremely low luminosity galaxies around 
the Milky Way. Another type of variable star, Cepheids, 
have already been used to id entify a candidate dwar f 
galaxy near the Galactic plane dChakrabarti ct ahll2015l ). 

In this paper, we investigate RR Lyrae stars in wide- 
held, time domain surveys as a possible tool for the dis¬ 
covery of dwarf galaxies with L > 1000 Lq. We are 
most interested in RRL as possible tracers of unseen 
dwarfs at low Galactic latitudes and/or with V-band cen¬ 
tral surface brightnesses fainter than 29 mag arcsec -2 . 
In Section [2J we review time domain surveys that are 
sensitive to RR Lyrae s tars at halo distances, and use 
the ELVIS simulations (jGarrison-Kinnncl et all 12014b ?) 
to quantify the fraction of predicted dark matter sub¬ 
halos within the grasp of planned surveys. In Section [3J 
we use the observed RR Lyrae populations of Milky Way 
dwarf galaxies to simulate reasonable representations of 
RR Lyrae in dwarf galaxies with a range of sizes and 


luminosities, apply a friends-of-friends algorithm to sim¬ 
ulated dwarfs to assess their detectability as groups of 
RR Lyrae stars, and apply a friends-of-friends algorithm 
to simulated stellar halos to assess the extent of inter¬ 
loper groups of held RR Lyrae stars. We summarize 
the detectability of our s imul ated dwarfs as groups of 
RR Lyrae stars in Section 13.51 and present the results of 
a friends-of-friends group hnder applied to existing RR 
Lyrae catalogs in Section [4] 

2. TIME DOMAIN SURVEYS SENSITIVE TO HALO RR 
LYRAE 

2.1. Summary of Survey Properties 

Numerous time domain surveys with published RR 
Lyrae catalogs are sensitive to RR Lyrae in the Milky 
Way’s halo. Existing RR Lyrae catalogs that cover 
more than 1% (^400 deg 2 ) of the sky at halo dis¬ 
tance s include: LINEAR (|Sesarl 120111: iPalaversa et al.1 
I20l3l) . PanSTARRS 1 - SPSS (lAbbas et all 120 lT ? 
La Silla QUEST (LSQ, IZinn et all I2014D . SEKBO 
(IKeller et all 120081) . the three component Catalina Sur¬ 
vey^ (|Drake et al.1 12009D . and Stripe 82 (ISesar et al.1 
I2010D . (Although SPSS Stripe 82 only covers 270 deg 2 
of sky, we include it in this list because of its excellent 
depth and RR Lyrae completeness.) 

Ongoing and future surveys that will yield catalogs of 
RR Lyrae stars at halo distances inc lude: the VISTA 
Variables in the Via Lactea (VVV, ISaito et all 120121 : 
IGran et~aLl I2015D . Palomar Transient Factory (PTF, 
and its successors Intermediate Palomar T ransient Fac¬ 
tory and Z wicky Transient Facility, lLaw et ahl 120091 : 
iBellml 12014 J. Surace private communication 2015), 
PanSTARRS 1 (B. Sesar private communicat ion 2015, N. 
Hern i tschek et al in preparation) , and LSST (llvezic et al.1 
120081 IVanderPlas fe Ivezidl2015l) . 

For reference, Table 1 summarizes the areal coverage, 
effective RR Lyrae depth (distance of ^50% complete- 

2 In cluding the Catalina Schmidt Survey (CSS, Drake et al.1 
2013a), the Sliding Spring Survey (SSS, Torrcalba ct a 120151) . and 
the Mount Lemmon Survey (MLS, [Drake ct al. 20 13b! ) 
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ness, in the absence of Galactic extinction), and maxi¬ 
mum RR Lyrae distance (when known) for each of these 
surveys. To calculate RR Lyrae distances from published 
apparent magnitudes, we assume [Fe/H] = -2.0 and My 
= 0.23([Fe/H]) + 0.93 which gives My ~ 0. 5 (iGhaboverj 
" 99!):. We also assume that My ~ M r , as in lSesar et al.l 
( 201311 . In many cases, the effective survey depth and 
maximum RR Lyrae distance are quite different. This 
is because the RR Lyrae completeness and purity drop 
precipitously roughly a magnitude brighter than most 
surv ey’s photometric li mits. However, in many cases 
(e.g. lDrake et al.l[2013al lbl). interesting Galactic structure 
work can be done at magnitudes fainter than the effective 
RR Lyrae depth. 

2.2. Fraction of Predicted Milky Way Sub-Halos 
Reached With RR Lyrae 

The time domain surveys that will be most relevant to 
the discovery of dwarf galaxies are those that are both 
wide-held and deep. To identify the time domain surveys 
with the greatest potential to yield new Milky Way dwarf 
galaxy discoveries, we quantify the fraction of predicted 
dark matter sub-halos that will be contained within RR 
Lyrae survey volumes as a function of areal coverage and 
depth. 

We use the publicPl ELVIS simulation suite 
(iGarrison-Kimmel et al.ll2014b ) to calculate the fraction 
of the dark matter sub-halos within 400 kpc that are 
enclosed within a given survey volume around the Milky 
Way, as a function of area and effective distance. ELVIS 
includes a set of 12 high-resolution, dark matter only 
simulations of Local Group-like environments simulated 
in a cosmological context. We made this prediction 
based on the average radial distribution of sub-halos 
around each of the Milky Way and M31 analogs in these 
12 Local Group-like simulations. We included sub-halos 
with Vpeak > 12 km s -1 and dhost < 400 kpc. 

Figure [2] shows this prediction, with contours outlin¬ 
ing the surfaces within which 5 - 80% of the predicted 
sub-halos are enclosed. The effective volumes of the sur¬ 
veys listed in Table 1 are also overplotted. As in §2.1, we 
assume My ~ M r ~ 0.5 for RR Lyrae and no Galactic 
extinction to go from effective V and r-depth to effec¬ 
tive distance. The estimated survey volumes in this sec¬ 
tion are therefore overestimates, because the effective RR 
Lyrae survey distances will be diminished in the Galactic 
plane. However, these estimated fractions of sub-halos 
to be included within each survey volume are sufficient 
to guide our exploration of this approach to dwarf find¬ 
ing. Of pre-LSST surveys, both PTF/iPTF/ZTF and 
PanSTARRS 1 have the grasp to map ^10-15% of sub¬ 
halos with RR Lyrae stars. 

Only LSST will be sensitive to RR Lyrae occupying a 
significant fraction of Milky Way sub-halos. LSST’s uni¬ 
versal wide-fast-deep (WFD) cadence will be sensitive 
to RR Lyrae well beyond the Milky Way’s viri al radius 
(|VanderPlas fe Ivezid 120151 : lOlusevi et all I20121) . includ¬ 
ing ~45% of predicted sub-halos. Because LSST is ex¬ 
pected to be much deeper than needed to reach the Milky 
Way’s virial radius, this 45% estimate will not be signif¬ 
icantly diminished by extinction. This is also a lower 
limit on the dwarf discovery space to be opened up by 

3 http://localgroup.ps.uci.edu/elvis/ 


LSST, because it doesn’t include (i) the volume between 
400 kpc < d < 600 kpc that will also be mapped by RR 
Lyrae with the WFD cadence, (ii) thousands of addi¬ 
tional square degrees to be observed with a non-universal 
cadence (e.g. the Northern Ecliptic Spur, Galactic Plane, 
and selected Deep Drilling Fields), and (iii) the fact that 
LSST’s Deep Drilling Fields could reveal RR Lyrae to d 
> 1.5 Mpc. 



Fig. 2.— The fraction of Milky Way sub-halos reached by RR 
Lyrae in current and planned time domain surveys, as a function 
of areal coverage and effective RR Lyrae detection distance. The 
contours show lines of sub-halo fraction contained within a given 
survey area and depth. These survey volumes are overestimates, 
because they do not account for Galactic extinction close to the 
plane. However, LSST is expected to be much deeper than needed 
to find RR Lyrae out to the Milky Way’s virial radius so this figures 
undersells the dwarf discovery volume to be enabled by LSST (at 
all but very low Galactic latitudes). 

This Figure does not consider any special model for 
how dwarf galaxies are predicted to occupy sub-halos. 
The fraction of Milky Way dwarf galaxies expected 
within these RR Lyrae survey volumes may, therefore, 
be somewhat different than the values show n in Figure [2l 
even i f we do reside in an LCDM universe. lHargis et all 
(1 2014 ?) used the ELVIS simulation suite to show that /> 
95% of the Milky Way’s dwarf galaxies are expected to 
reside at d > 50 kpc, for several reasonable models used 
to populate dark matter sub-halos with dwarf galaxies. 
However, the fraction of MW dwarf galaxies expected to 
reside at d > 100 kpc ranges from 50% - 80%, depending 
on model. 

This Figure also doesn’t consider the detectability of 
any of these sub-halos (dwarf galaxies) as distinct groups 
of RR Lyrae stars. In Section 13.41 we will show that de¬ 
tecting dwarf galaxies as groups of RR Lyrae stars that 
are distinct from the field becomes easier at larger dis¬ 
tances, where few field RR Lyrae are expected. There¬ 
fore, deeper surveys for RR Lyrae will generally be more 
fruitful for dwarf galaxy discovery than wider surveys. 

3. RR LYRAE AS A TOOL FOR DWARF GALAXY 
DISCOVERY 

There are two things that affect a dwarf galaxy’s de¬ 
tectability as a group of RR Lyrae: (i) whether a given 
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TABLE 1 

Sensitivity of Time Domain Surveys to RR Lyrae Stars 


Survey 

Coverage 

deg 2 

effective 3, 
mag limit 

d e// a 
kpc 

<i b 

'-*max 

kpc 

Public RRL 
Catalog 

Completeness to RRab 
and Comments 

LINEAR 

10,000 

r~17.5 

~25 

32 

Y 

~80% for r<17.2 

SDSS-PanSTARRS-Catalina 

14,000 

g~17.8 

~28 

37 

Y 

~50% for V<17.8 

CRTS CSS 

~20,000 

V~18.0 

~30 

80 

Y 

~50% at V = 18.0 

CRTS SSS 

14,800 

V~18.0 

~30 

60 

Y 

~60% at V = 18.0 

CRTS MLS 

~1,000 

V~20.0 

~75 

125 

Y 

~50% at V = 20.0 

SEKBO 

1675 

V~18.5 

~40 

95 

Y 

~60% for V<18.5 

La Silla QUEST 

840 c 

r~19.0 

~50 

95 

Y 

>70% for V <19. 

Stripe 82 

270 

r~21.0 

~120 

150 

Y 

99% to 120 kpc 

vvv 

562 

Ks~17 

~14 

- 

N 

ongoing 

PTF/iPTF/ZTF 

30,000 d 

r~20.0 

~80 

- 

N 

ongoing 

PanSTARRS 1 

30,000 

r~20.5 

~100 

- 

soon 

data release in preparation 

LSST Wide-Fast-Deep e 

~18,000 

r~24.5 

~600 


TBD 

future, 100% complete at r~24.5 
images and photometry will be public 


Note. —References: LINEAR -ISesar 

(2011 

): Palaversa et al.| (20131); PanSTARRS 1 - SDSS - 

Abbas et al. (2014b; CRTS CSS/MLS 

- Drake et al. (20091), iDrake et al. (2013a 

), and 

Drake et al. (2013bl): CRTS SSS -iTorrealba et al. 


201E 

); SEKBO - Keller et al.l (2008|): 

LSQ - Zinn et al. (20141): Stripe 82 - Sesar et al.I (120101); VVV - iSaito et al. (20121); [Gran et al. 

201.E 

); PTF/iPTF/ZTF - Law et al. 


<20031 . 155153 <20lib J . Su race (private comm u nication, 2015); Pan STARRS 1 - B. Sesar (private communication 2015), Hernitschek et 
al. (in preparation); LSST- [lvezic et al ] d2008l) , [VanderPlas Ivezid J2015I) 


3 An effective magnitude and distance limit for RR Lyrae detection, based on the magnitude at which the RR Lyrae sample is reasonably 
complete, assuming no extinction and My = 0.5. “Reasonably” ranges from ~50% to ~90% complete, depending on the information 
provided in the cited paper(s). 

k The maximum distance of a cataloged RR Lyrae, for surveys with publicly available catalogs. 

c 840 deg 2 is the “Region I” area discussed in IZinn et al . <20Tl . The LSQ will eventually survey 15,000 deg 2 , 

http: //hep.yale.edu/lasillaquest 

The upper limit of the survey area to be covered by iPTF (J. Surace, private communication 2015) and/or ZTF. 
e We specify “Wide-Fast-Deep” here to distinguish the part of the LSST survey that will be observed with a universal cadence, as 
described in llvezic et al.l d2008l) , from other survey extensions such as Deep Drilling Fields and the Northern Ecliptic Spur. This is an 
underestimate of the area over which LSST data will be able to recover RR Lyrae. 


dwarf galaxy has enough RR Lyrae stars with a compact 
enough spatial distribution to be identified as a group, 
and (ii) whether the number of interloper field RR Lyrae 
groups is low enough such that true dwarf galaxy groups 
are distinct. To determine (i), we assume a population of 
RRab stars similar to those in known Milky Way dwarf 
galaxies (Section [34]) and simulate the spatial distribu¬ 
tion of RRab stars in dwarf galaxies with -10 < My 
<-2.5 and 10 pc < rhaif < 1000 pc (Section 13.21) . To 
determine (ii), we apply a 2D friends-of-friends (FOF) 
group-finding algorithm (Section 13.01) to simulated RR 
Lyrae catalogs to study the expected demographi cs o f 
interloper “groups” of field RR Lyrae stars (Section l3~il) . 
Informed by our investigation of the field RR Lyrae pop¬ 
ulation, we apply an FOF algorithm with a linking length 
of 500 pc to the simulated dwarf galaxies to measure their 
detectability as a function of linking length, luminosity, 
and half-light size (Section 13.51) . 

3.1. Known RR Lyrae Populations of MW Dwarf 
Galaxies 

The Milky Way’s dwarf spheroidal and ultra-faint 
dwarf satellites span four orders of magnitude in lumi¬ 
nosity, including dwarfs with only several hundred So¬ 
lar luminosities. Despite the low luminosities of many 
Milky Way dwarfs, at least one RRab star has been 
found in every dwarf for which ther e is p ubli shed time- 
seri es observations. For e xample, ISirnon et al.1 (|20111 ) 
and IBoettcher et al.1 (|2013f) discovered one RR Lyrae in 
each of Segue 1 and Segue 2 respectively, which are 
the two faintest M W dwarf galaxies currently known 
(|McConnachid[20l^ ) . On the other end of the MW dwarf 
spectrum, more than 100 RR Lyrae have been discovered 
in each of Draco, Leo II and Sculptor dBonanos et al.l 



Fig. 3.— The observed number of RRab stars in Milky Way dwarf 
galaxies as a function of V-band absolute magnitude. Three lines 
of constant specific frequency, S'rr = -/Vrr * io°- 4 ( 7,5 + m v) 5 are 
overplotted. The specific frequencies of the most luminous dwarfs 
on this figure are lower limits, because of incompleteness. Specific 
frequencies of 20 - 100 are the most representative of dwarf galaxies 
in the magnitude range we are most interested to find. 

12004ISiegel fc Maiewskill20t)0HKahiznv et al.lll995t) . 

Table m summarizes the known number of RRab and 
RRc stars in 17 Milky Way dwarf galaxies and Figure [3] 
shows the observed number of RRab stars as a function of 
each Milky Way dwarf satellite’s absolute V-band mag¬ 
nitude. We focus Figure [3] on only RRab stars for three 
reasons: RRab are the most numerous type of RR Lyrae 
star, RRab stars are easier to identify than RRc stars 
because of their distinctly shaped light curves and larger 
amplitudes, and RRab stars have well calibrated [Fe/H]- 
luminosity relationships. 

We overplot several v alues of specific frequency , using 
the parameterization of iMackev fc Gilmore! (120031) : 
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TABLE 2 

Observed Milky Way Dwarf RR Lyrae Populations 


MW Dwarf 

N RRab 

N RR C 

RR c /RR a 6+c 

Reference 

Fornax 

396 

119 

0.23 

Bersier & Wood (2002) 

Sculptor 

134 

88 

0.40 

Kaluznv et al. (1995) 

Draco 

214 

30 

0.12 

Kinemuchi et al. ?2008) 

Leo II 

106 

33 

0.24 

Siegel & Maiewski (2000) 

Ursa Minor 

47 

35 

0.43 

Nemec et al""£l988) 

Leo I 

47 

7 

0.13 

Held et al. (2000^ 

Sextans 

26 

10 

0.27 

Mateo et al. (1995) 

Canes Venatici I 

18 

5 

0.22 

Kuehn et al. (2008) 

Bootes I 

7 

8 

0.53 

Siegel (20011 

Hercules 

6 

3 

0.33 

Musella et al. (2012) 

Ursa Major I 

5 

2 

0.29 

Garofalo et al. (2013) 

Leo IV 

3 

0 

0.0 

Moretti et al. (20091 

Coma Berenices 

1 

1 

0.5 

Musella et al. (2009) 

Canes Venatici II 

1 

1 

0.5 

Greco et al. (20087 

Ursa Major II 

1 

0 

0.0 

DairOra et al. (2012) 

Segue 2 

1 

0 

0.0 

Boettcher et al. ^2(^5) 

Bootes II 

1 

0 

0.0 

Sesar et al. (2014) 

Segue 1 

1 

0 

0.0 

Simon et alM25ff) 


Srr = Nrr * io°- 4 ( 7 - 5 + M v) ( 2 ) 

where Nrr is the number of RR Lyrae in the dwarf 
galaxy and specific frequency, Srr, is the number of 
RR Lyrae stars observed per 10,000 L 0 . This Figure 
highlights a factor of 20 difference in observed specific 
frequency between the least luminous satellites (Srr ~ 
100) and the more luminous satellites (Srr ~ 5). In¬ 
completeness in the RR Lyrae census in the brighter 
dwarf galaxies likely accounts for some of this difference. 
Differences in stellar populations may account for the 
rest of the difference. For example, there is an observed 
stellar mass-metallicity relation and an observed stellar 
mass-SFH relation (lower mass dwarfs seem to stop form¬ 
ing stars earlier iKirbv et al.11201.11 iBrown et all 12014 
IWeisz et al.ll2015l) . 

We adopt an Srr of 20 - 100 to be representative of the 
low luminosity dwarf galaxies that RR Lyrae may help 
discover. This Srr gives a lower limit to the number 
of RR Lyrae that would be cataloged by studies that 
are also sensitive to RRc stars. Known RRc fractions in 
the Milky Way’s dwarfs range from ~0.2 - 0.5. The total 
number of RR Lyrae hosted by a dwarf galaxy is therefore 
30% - 100% larger than expected for RRab stars alone. 

3.2. Simulating the RR Lyrae Populations of Dwarfs 

We simulated the RR Lyrae spatial distributions of 
13,200 dwarf galaxies with -10 < My < -2.5 and 10 pc 
< half-light radius <1000 pc. Our idealized sim ulati ons 
assume 100% RR Lyrae completeness. In Section Ebll we 
describe our simplistic RR Lyrae group finding approach. 
It is based on a physical, rather than an angular, linking 
length so a dwarf’s distance doesn’t affect how detectable 
it is as a group of RR Lyrae. (Although more distant 
dwarfs will be easier to pick out from the field RR Lyrae 
population - see Section 13.41) . 

13,200 dwarfs were simulated for each of two assumed 
RR Lyrae specific frequencies, Srr = 20 and 100. The 
number of RR Lyrae assigned to each simulated dwarf 
was drawn from a Poisson distribution with an average 
given by Srr/ 10^ 0 4 ^ 7 5+Mv ^. The simulated stars were 
spatially distributed using an exponential distribution. 
These simulations only account for the RRab population 


of dwarf galaxies, so underestimate the total number of 
RR Lyrae stars. 

3.3. Friends-of-Friends Approach to RR Lyrae Group 

Finding 

We apply a simplistic 2D FOF algorithm (IDL’s 
spheregrou]Q) for RR Lyrae group finding to each sim¬ 
ulated dwarf, as well as to simulated MW stellar halo 
models (Section l3~4ll and public RR Lyrae catalogs (Sec¬ 
tion SI). We apply the algorithm on stars in overlapping 
bins of one apparent magnitude. When initially inves¬ 
tigating this method for detecting dwarfs, we tested 2D 
linking lengths between 50 pc and 500 pc (the physi¬ 
cal size range observed for MW dwarf galaxies) to inves¬ 
tigate the trade-off between increased recoverability of 
low surface brightness dwarfs as RR Lyrae groups and 
increased contamination by “groups” of field RR Lyrae 
stars when using larger linking lengths. We ultimately 
found that even a 500 pc linking length results in an ac¬ 
ceptable level of field RR Lyrae group contamination at 
outer halo distances (Section [33}, so we will focus on the 
500 pc linking length results for the rest of this paper. 

3.4. The Predicted Population of Interloper RR Lyrae 

Groups 

The extent to which spatial associations of unbound 
field RR Lyrae stars will contaminate RR Lyrae-based 
searches for bound MW dwarf galaxies will be a function 
of distance, linking length, and the number of RR Lyrae 
required to define a group. We focus our investigation 
of interloper RR Lyrae groups at d > 50 kpc, both be¬ 
cause there are a large number of field RR Lyrae with d 
< 50 kpc and because ^95% of the MW’s dark matter 
sub-halos are predicted to have d > 50 kpc (Figure [2]). 
Although we tried linking lengths of 50 pc - 500 pc in our 
analysis, we only present the number of interloper field 
groups identified with a 500 pc linking length because, as 
we show below, that linking length yields a small number 
of interloper field RR Lyrae groups. 

There is currently no observed sample of RR Lyrae 
stars that is both wide-held and deep enough to thor¬ 
oughly investigate the demographics of held RR Lyrae 

4 http://spectro.princeton.edu/idlutils_doc.html 
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Fig. 4.— The spatial distributions of RR Lyrae analog stars in 
semi-analyti c stellar halo models applied to the Aqu arius A - D 
simulations (Springel ct al. 2008; Lowing et al. 2015). Stars were 
selected by applying criteria of (i) 6100 K < T e ff < 7400 K, (ii) 
2.5 < log(g) < 3.0, (iii) unbound to any galaxy, and (iv) more 
than 50 kpc from the “S un” to th e simulated stellar halo catalogs 
generated by I Lowing et all (120151) . The highly structured spatial 
distributions, and total number of RR Lyrae stars, are both in 
agreement with past predictions and current observations. 


at cl > 50 kpc (although we use the Stripe 82 catalog as 
a sanity check below). We therefore turn to simulations 
to describe the expected properties of RR Lyrae in the 
MW’s stellar halo. (See Section 5 for caveats related to 
possible thi ck d isk conta mination in the Galactic plane.) 

iLowing et al.l (|2015f) applied th e semi-analytic 
galaxy formation model, GALFORM dGole et, all 120001 
iFont et ahl 1201 lh , to five dark matter only simulations 
of Milky Way analogs. They used the Aquarius A - 
Aq uarius E halos fro m the Aquarius simulation project 
(jSpringel et al.l 120081 1. These are all relatively is olated 
halos with M to t ~ 10 12 M 0 . ILowing et al.l (12015H 
combi ned the particle tagging approach of iCooper et al l 
rt2010h with stellar popula tion synthesis and p hase space 
sampling (ENBID, iSharma fc Steinmetzl 1200(1) to gen¬ 
erate full catalogs of halo stars in each simulation. Like 
any simulation, this method has some limitations. For 
example, it only accounts for the accreted population 
of halo stars, it doesn’t include the lowest luminosity 
dwarf galaxy accretions (which contribute little to the 
halo), and stars are initially given the phase space 
distribution of the dwarf galaxy dark m atter parti cles 
onto which they are tagged. However, ILowing et all 
(12015T) demonstrated the efficacy of these stellar halo 
catalogs to provide faithful realizations of the majority 
of the accreted stellar halo at d > 20 kpc. Moreover, we 
checked that the simulated z = 0 dwarf galaxy satellites 
have specific frequencies of RR Lyrae analog stars (as 
selected below) that are similar to those in observed 
dwarf galaxies. For the Aquarius B and D simulations, 
the Srr of well resolved dwarfs with 50 kpc < d < 300 
kpc ranged from ~10 to 100, with most between 20 - 70 
(a bit higher than observed for the Milky Way’s more 
luminous dwarfs, as described in Section 3.1). 


To select RR Lyrae analog stars from these catalog^, 
we applied selection criteria of (i) 6100 K < T e g < 7400 
K, (ii) 2.5 < log(g) < 3.0, (iii) unbound t o any galaxy , 
and (iv) more than 50 kpc from the “Sun” (lSmithlll995li . 

The number of analog RR Lyrae stars satisfying these 
criteria ranges from ~1,500 - 4,000 (Aquarius B and E) 
to ~ 35,000 - 42,000 (Aquarius A, C, and D). These num¬ 
bers are consistent with the number observed in Stripe 
82, the only public RR Lyrae catalog that is complete 
to d ~ 100 kpc. There are 55 Stripe 82 RR Lyrae with 
50 kpc < d < 100 kpc. A crude area correction yields 
an expected ~8,100 RRL in the MW with 50 kpc < d < 

100 kpc, right in the middle of the predicted number in 
Aquarius A-E (622 to 39,738). 

The spatial distributions of these simulated RR Lyrae 
field populations are also qualitatively similar to previ¬ 
ous N-body+semi-analytic predictions of t he stel lar halos 
of MW-mass galaxies (e.g. I.Tohnston et aDl2008 L Figure 
0] shows the spatial distributions of Aquarius A-D, in 
Cartesian coordinates, with the x-axis aligned with the 
major axis of the simulated dark matter halo. (We ex¬ 
clude Aquarius E because it has a very small number of 
RR Lyrae with d > 50 kpc and is qualitatively consistent 
with Aquarius B.) These figures show that the predicted 
RRL distributions beyond 50 kpc are highly structured, 
with significant differences between the halos, and few 
RR Lyrae with d > 100 kpc (<; one per 10 deg 2 on aver¬ 
age, 965 - 5,739 total). 

We applied the FOF algorithm with a 500 pc linking 
length, as described in Section 13.21 to each of the five 
simulated stellar halos. Depending on whether we im¬ 
plemented a group threshold of 2 or 3 RR Lyrae stars, 

Nrr, each simulation contained dozens to a couple thou¬ 
sand field RR Lyrae “groups” with d > 50 kpc. Although 
this sounds like a large number of field RR Lyrae groups, 
these groups are highly spatially clustered around a small 
number of rich stellar streams, with few isolated field 
RR Lyrae groups. The number of field RR Lyrae groups 
found in the simulations with 50 kpc < d < 100 kpc 
is consistent with the number found in the 280 deg 2 of 
Stripe 82. Five groups of 2 stars are found in Stripe 82 
and no groups were found with three or more stars. (See 
Section |4] for more discussion). Assuming Stripe 82 is 
representative of the typical field halo RR Lyrae popula¬ 
tion around the Milly Way, the total number of groups 
expected is ~700. This is in between the small number 
found in Aquarius B and E (~200 - 300) and the number 
found in Aquarius A, C, and D (~2000 - 3000). 

We first focus on the predicted number and spatial dis¬ 
tribution of field RR Lyrae groups with d > 100 kpc in 
the simulations, which includes the vast majority of the 
MW’s virial volume and ~80% of predicted dark mat¬ 
ter sub-halos. Figure [5] shows the distribution of such 
RR Lyrae analog stars in Aquarius B and Aquarius D 
in grey. Groups identified with two or more RR Lyrae 
stars are overplotted in red. At these large distances, 
there are fewer than 5 field RR Lyrae groups are isolated 
from obvious large-scale halo overdensities. Dwarfs at d 
> 100 kpc can therefore be identified with little contam¬ 
ination (<; 1 per 8000 deg 2 on average over large areas) 
as groups of 2 or more RR Lyrae linked within 500 pc, 

5 available at http: //galaxy-catalogue. dur. ac .uk: 8080/StellarHalo/ 








































Dwarf Galaxy Searches with RRL 


7 


by ignoring RR Lyrae groups embedded near large-scale 
halo overdensities (which appear to occupy <1/6 of the 

sky)- 

We next focus on the predicted number and spatial 
distribution of field RR Lyrae groups with 50 kpc < d < 
100 kpc, which includes a few percent of the MW’s virial 
volume and ~15% of predicted dark matter sub-halos. 
Figure [6] is the same as Figure [5] but for 50 kpc < d < 
100 kpc and shows that many more field RR Lyrae stars 
exist with 50 kpc < d < 100 kpc, with large-scale halo 
overdensities occupying up to 1/2 of the sky in projec¬ 
tion. The top two panels of Figure [5] show the distri¬ 
bution of groups with two or more RR Lyrae and the 
bottom two panels show the distribution of groups with 
three or more RR Lyrae. There are as many as a few 
dozen relatively isolated groups found with two or more 
RR Lyrae, but less than 10 relatively isolated groups of 
three or more RR Lyrae. With fewer than 10 relatively 
isolated field RR Lyrae groups, dwarfs at 50 kpc < d < 
100 kpc could be identified with little contamination (<; 

1 per 4000 deg 2 on average over large areas) as groups 
of 3 or more RR Lyrae linked within 500 pc, by ignoring 
RR Lyrae groups embedded near large-scale halo over¬ 
densities (which occupy <1/2 of the sky). Groups with 
only 2 stars can also reveal new dwarf galaxies with 50 
kpc < d < 100 kpc when quite isolated from other field 
RR Lyrae clusterings. 

3.5. The Detectability of Milky Way Dwarfs Using RR 

Lyrae 

We focus our discussion of detectability on dwarfs that 
would be identified as groups of two or more RR L yrae 
stars, with a 2D linking length of 500 pc. In Section I7TTT1 
we demonstrated that we expect few isolated, field RR 
Lyrae groups of two of more RR Lyrae in the outer halo 
(d > 100 kpc). These distances include the vast majority 
of the MW’s halo volume and predicted dwarf galaxies. 
Groups of two or more RR Lyrae stars at nearer distances 
(50 kpc < d < 100 kpc) that are quite isolated from visi¬ 
ble RR Lyrae structures will also distinctly identify dwarf 
galaxies. In the discussion below, the luminosity detec¬ 
tion limits would be ^0.5 mags brighter with a threshold 
of three RR Lyrae, as might be desired at these nearer 
distances (50 kpc < d < 100 kpc) to reduce the num¬ 
ber of expected field RR Lyrae groups that are spatially 
coincident with diffuse halo structures. 

Figure [7] shows the detectability of simulated dwarfs 
as groups of two or more RR Lyrae stars, as a func¬ 
tion of size, absolute magnitude, and RR Lyrae specific 
frequency, Srr. The white line traces the sizes and lumi¬ 
nosities of dwarfs that are identified as RR Lyrae groups 
50% of the time. The left (right) panel shows the de¬ 
tectability of dwarfs with Srr, = 20 (100). These specific 
frequencies are similar to those observed in dwarfs more 
(less) luminous than My ~ -4.5 mag (see Section [3711) . In 
fact, all dwarfs less luminous than My ~ -6 have RRab- 
only specific frequencies of at least 20. The simulations 
informing the left panel therefore include the fewest RR 
Lyrae stars expected in ultra-faint dwarfs, yielding con¬ 
servative detection limits at those luminosities. 

For a specific frequency of 20, dwarfs similar to those 
known, with My brighter than ~ -5 mag, will be iden¬ 
tified as groups of 2 or more RR Lyrae stars. Assuming 


Srr = 20, low surface brightness dwarfs (at least as faint 
as Hv,o ~31 mag arcsec -2 ) that are more luminous than 
My = -6 mag will also be detectable as groups of two 
or more RR Lyrae stars. Even with our conservative 
estimate of the number of RR Lyrae to belong to these 
dwarfs, this demonstrates that groups of two or more RR 
Lyrae will distinctly reveal (i) dwarfs brighter than My 
= -5 mag that are similar to those known, but at low 
Galactic latitudes and (ii) slightly more luminous dwarfs 
at previously unexplored surface brightnesses. 

The detectability of dwarfs with luminosities fainter 
than My = -5 mag will depend sensitively on their RR 
Lyrae specific frequency. Assuming a specific frequency 
of 100 for the least luminous dwarfs (e.g. My less lumi¬ 
nous than -4.5 mag), the right panel demonstrates that 
dwarfs similar to those known and with luminosities as 
low as My = -3.2 mag will be detected as groups of two 
or more RR Lyrae stars. Dwarfs will also be detected at 
these low luminosities, with surface brightnesses as faint 
as 31 mag arcsec -2 . This demonstrates that groups of 
two or more RR Lyrae will also distinctly reveal dwarfs 
with -4.5 mag < My < -3.2 mag that are (i) at low Galac¬ 
tic latitudes and/or (ii) at previously unexplored surface 
brightnesses (depending on RR Lyrae specific frequency). 

Although RR Lyrae will open up significant new dwarf 
discovery space, there are some dwarfs that RR Lyrae 
will not efficiently detect. These include those less lu¬ 
minous than My = -3.2 mag unless they have an RR 
Lyrae specific frequency higher than yet observed, dwarfs 
fainter than My = -5 mag with specific frequencies less 
than 20, and dwarfs with -6 mag < My < -5 mag and 
lower surface brightness than 31 mag arcsec -2 . 

4. SEARCHING FOR RR LYRAE GROUPS IN PUBLIC 
SURVEY DATA 

Although most published public RR Lyrae catalogs are 
severely incomplete for d > 50 kpc (and all but Stripe 82 
are severely incomplete for d > 100 kpc), several of them 
include RR Lyrae with d > 80 kpc (CRTS - the CSS 
and MLS, SEKBO, LSQ, and Stripe 82 - see Table 1 and 
references therein). It is worth searching these catalogs 
for isolated groups of RR Lyrae, because very low sur¬ 
face brightness or low Galactic latitude dwarfs with the 
luminosity of (for example) Sextans or Draco, could be 
identified as an RR Lyrae group even with 10% RR Lyrae 
completeness. We therefore apply the FOF group fi nder 
with a 500 pc linking length, as described in Section ITTT7T1 
to these catalogs. As before, we only consider RR Lyrae 
with d > 50 kpc and, when possible, we included only 
cataloged RRab stars. To consistently estimate distances 
to the RR Lyrae in each of these su rveys, we assumed My 

M r = 0.5 mag. Wc also applied ISchlearel et, all (|1998 1 
corrections for dust extinction to the LSQ and CRTS 
magnitudes, and applied the correction for dust extinc¬ 
tion as provided in the SEKBO catalog. The Stripe 82 
RR Lyrae catalog provided with an extinction correction 
applied. 

Figure [5] shows the spatial distributions of RR Lyrae 
and RR Lyrae groups in the CRTS, LSQ, SEKBO, and 
Stripe 82 catalogs. Grey symbols show the spatial distri¬ 
bution of all RR Lyrae in these catalogs, black symbols 
show the RR Lyrae contained within each of three dis¬ 
tance bins (50 kpc < d < 60 kpc, 60 kpc < d < 70 kpc, 
and d > 70 kpc), and colored symbols show groups of 
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Fig. 5.— The spa tial distribution of RR Lyrae analog star s at d > 100 kpc in semi-analytic stellar halo models applied to the Aquarius B 
and D simulations IlSpringel et aljf2008t ILowing et al.ll20 1 5'). shown in Equatorial coordinates. Red symbols show friends-of-friends groups 
of two or more simulated field halo stars, identified with a linking length of 500 pc. At these outer halo distances, very few isolated groups 
of field RR Lyrae stars are expected to contaminate searches for dwarf galaxies as groups of RR Lyrae stars. 



Fig. 6.— Top panels: The spatial distributio n of RR Lyrae analog stars at 50 kpc < d < 100 kpc in semi-analytic stellar halo models 
applied to the Aquarius B and D simulations IlSpringel et all 120081 ; ILowing et al .1120 15 ). shown in Equatorial coordinates. Red symbols 
show groups of two or more simulated field halo stars, identified with a linking length of 500 pc. Although there are as many as a few 
dozen relatively isolated field groups or two or more RR Lyrae in this distance range, groups with only 2 stars will distinctly reveal new 
dwarf galaxies with 50 kpc < d < 100 kpc when quite isolated from other field RR Lyrae spatial clusterings. Bottom panels: Same as top 
panels, but for groups of three or more RR Lyrae stars. Fewer than 10 relatively isolated groups of field RR Lyrae are found with this 
higher threshold. Dwarfs in this distance range could therefore be identified as groups of three or more RR, Lyrae stars with < 1 interloper 

group per 4000 deg 2 , on average. 


two or more RR Lyrae. 

A visual examination of the RR Lyrae distributions in 
Figure |8] reveals several known halo overdensities. In the 
top and middle panels, the Sagittarius (Sgr) leading arm 
is clearly seen with 180° <; RA <, 240°. By compari¬ 
son with iTorrealba et all (|2015l ) the diffuse overdensity 
with 150° RA < 180° and 0° Dec < 30° in the 
top and middle panels is also likely Sgr stream mate- 
rial. The Gemini stream (a likely distan t wrap of Sgr, 
iDrake et al.ll2013bl : iBelokurov et al.l|2014f l is clearly visi¬ 
ble in the bottom panel of Figure[5]at 100° < RA <120° 
and Dec ~ 20°, with a possible extension to the stream¬ 


like feature seen with 180° <; RA <; 210° and slightly 
negative declination. The Pisces overdensity is also seen 
in the bottom panel of F i gure p near RA ~ 350° an d 
Dec ~ 0°. SeelSesar et all (l20071blWatkins et all (I2009T) . 
and lSesar et all (120101) for the discovery of and follow-up 
discussions about Pisces. 

4.1. Relatively Isolated RR Lyrae Groups in Public 
Catalogs 

We selected 18 groups of two or more RR Lyrae that 
are not embedded in high surface density halo structures, 
for further inspection. We didn’t require “relatively iso- 
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Fig. 7. — The detectability of simulated dwarfs as groups of 2 or more RR Lyrae stars identified with a 500 pc linking length, as a 
function of size, My, and RR Lyrae specific frequency, Srr. A group threshold of 2 is most appropriate for d > 100 kpc, or for 50 kpc < 
d < 100 kpc if some contamination in the sample of dwarf candidates is acceptable or if groups found near large overdensities are excluded 
(see Figure Blue circles show known (filled points) and candidate (open points) MW dwarfs. The white line traces the sizes and 
luminosities of dwarfs that are detected as RR Lyrae groups 50% of the time. The results in the left (right) panel assumes an RRab star 
specific freq uency of 20 (100). Dwarfs with relatively lower specific frequencies are likely those more luminous than My = -4.5 mag (see 
Section CD}' Identifying dwarfs as groups of 2 or more RR Lyrae stars will open up significant discovery space at low Galactic latitudes 
and low surface brightnesses. 


TABLE 3 

Relatively isolated RR Lyrae groups in public catalogs 


RA 

Dec 

Avg M v 

Distance 3, 

Nrr 

Size (pc) b 

Survey 

Classification 

153.17955 

-2.04230 

20.2 

84 

2 

328 

LSQ 

Sextans 

209.85728 

14.40900 

19.4 

59 

4 

674 

CRTS 

Bootes I 

120.40614 

12.16494 

19.3 

57 

2 

150 

SEKBO 

likely Sgr 

140.01188 

22.67615 

19.3 

57 

2 

513 

CRTS 

likely Sgr 

168.27063 

10.05624 

19.1 

52 

2 

505 

CRTS 

likely Sgr 

178.85546 

13.87637 

19.4 

59 

2 

107 

CRTS 

likely Sgr 

180.65865 

23.47226 

19.4 

58 

2 

406 

CRTS 

possible Sgr? 

218.35120 

-14.26636 

19.5 

61 

2 

172 

SEKBO 

possible Sgr? 

241.24207 

-23.50549 

19.3 

55 

3 

638 

SEKBO 

likely Sgr 

316.12666 

1.11801 

19.2 

53 

2 

382 

CRTS 

possible Sgr? 

331.79452 

0.24400 

19.1 

51 

2 

527 

Stripe 82 

possible Pisces/Sgr? 

344.55949 

0.27140 

19.6 

64 

2 

523 

Stripe 82 

possible Pisces/Sgr? 

345.44403 

-9.98595 

19.3 

55 

2 

415 

SEKBO 

possible Pisces/Sgr? 

148.68102 

16.03975 

19.9 

76 

2 

257 

CRTS 

likely Gemini 

150.01648 

22.52547 

19.9 

74 

2 

360 

CRTS 

possible Gemini 

200.89584 

-5.36901 

20.1 

83 

2 

553 

SEKBO 

possible Gemini 

266.29432 

36.87977 

19.1 

52 

2 

531 

CRTS 

unknown 

354.56833 

25.72140 

19.3 

57 

2 

165 

CRTS 

unknown 


Note. — Groups of 2 of more RR Lyrae stars identified with a 500 pc linking length, with d > 50 kpc, 
and not spatially coincident with dense RR Lyrae structures. 

a If an RRL group fell right on the edge of a distance bin, its recorded distance here may differ slightly from 
its location in Figure [si 

k Sizes of groups of 2 RR Lyrae may exceed the linking length (500 pc) by a small amount due to our method 
of group detection: When converting angular separations to physical separations when group finding, we use 
the middle of the distance bin being searched. When calculating the physical size of each group, we utilize 
the mean distance of stars belonging to the group. 

lated” to be all that isolated when selecting these groups, 
in order to be fully inclusive of possibly interesting groups 
while knowingly including interloper field halo groups. 

While 18 RR Lyrae groups is loosely consistent with ex¬ 
pectations for the number of field RR Lyrae groups found 
in Section lTdl the observed and predicted numbers aren’t 
exactly comparable because of the significant incomplete¬ 
ness at d > 50 kpc in these RR Lyrae catalogs. 

The properties of these groups are presented in Ta¬ 
ble [3] and their positions are highlighted in Figure [8] as 
the colored stars with boxes around them. Two of these 
18 groups are associated with known MW dwarf galax¬ 
ies: the Bootes I dwarf and Sextans. These results are as 


expected, because no Milky Way dwarf with more than 1 
known RR Lyrae star resides within any public RR Lyrae 
survey’s effective survey volume. Although RR Lyrae are 
cataloged beyond these effective distances, completeness 
rapidly decreases and may be variable with location. Sex¬ 
tans was found as a group of 2 RRab stars and resides 
with d = 86 kpc in the MLS footprint, a survey with 
a 50% completeness distance of 75 kpc and a maximum 
detection distance of 125 kpc. Bootes I was found as a 
group of 4 RRab stars and resides at 66 kpc in the CSS 
footprint, a survey with a 50% completeness distance of 
30 kpc and a maximum detection distance of 80 kpc. It 
is therefore mildly surprising that 4 of its 7 of its RRab 
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Fig. 8. — Spatial distribution of RR Lyrae in the CRTS, LSQ, 
SEKBO, and Stripe 82 public catalogs, in Equatorial coordinates. 
The grey symbols in all three panels demonstrate the combined 
survey footprint, with the spatial distribution of all RR Lyrae in 
these surveys. Each panel also shows RR Lyrae (black symbols) 
and groups or two or more RR Lyrae (colored symbols) in different 
distance slices. Groups liberally classified as ’’relatively isolated” 
have boxes around them; 14 of of these likely belong to known halo 
structures, two to the Bootes I and Sextans dwarf, and two are 
unknown. Top panel: 50 kpc < d < 60 kpc; Middle panel: 60 kpc 
< d < 70 kpc; Bottom panel: d > 70 kpc). To consistently estimate 
distances to the RR Lyrae in all of these surveys, we assumed My 
~ M r = 0.5 mag. There are thus modest differences between the 
distances shown here and shown in the RR Lyrae survey papers. 

stars were detected as a group. Ursa Minor (77 kpc) and 
Draco (76 kpc) are the only other dwarfs known with 
more than 2 RR Lyrae stars that lie within the max¬ 
imum detection distance of a public RR Lyrae survey 
(CSS), but they are only a couple of kpc within that 
maximum distance and effective RR Lyrae completeness 
at the locations of those two dwarfs could easily be 0% 
at those distances. 

14 of these 18 groups appear to be associated with 
known halo structures. Nine of the highlighted groups 
with 50 kpc < d < 60 kpc are likely t o be associated 
with Sgr debris. We use Figure 16 of iTorrealba et all 
(|2015f) as a guide to classify the four with 120° <; RA 
<; 180° and Dec ;> 0° as likely Sgr debris, the one at 


RA ~ 240° and Dec ~ -25° as likely Sgr debris., and 
the one at RA ~ 180° and Dec ~ 25° as possible Sgr 
debris. The three groups with 315° <; RA < 330° and 
Dec ~ 0° are more ambiguous, and may belong to bits 
of Sgr or Pisces (for example, see the Pisces extension 
visible in Figure 10 of lSesar et al.ll2010h . Similarly both 
highlighted groups between 60 kpc < d < 70 kpc, could 
potentially be associated with Sgr debris, though the one 
with R A ~ 345° could also be part of the Pisces overden¬ 
sity. Finally, all of the isolated groups at d > 70 kpc (not 
including the Sextans dwarf) may be associated with the 
Gemini structure. 

The two structures we classify as “unknown” reside at 
RA ~ 355° and Dec ~ 26° and RA ~ 240° and Dec ~ 
40°. Although the former lies in the general direction of 
the Tri-And overden sity, it is 20 kpc more distant (e.g. 
iSheffield et al.ll20l4lh Because these two groups are at 
d < 60 kpc and reside close to other halo RR Lyrae, it 
is likely that they are field RR Lyrae stars. However, it 
would be worth following-up those regions of sky with 
imaging and/or a more complete RR Lyrae catalog. 

5. DISCUSSION AND CONCLUSIONS 

In this paper, we provide a proof-of-concept that 
groups of two or more RR Lyrae stars can distinctly re¬ 
veal MW dwarf galaxies at d > 50 kpc, and at Galac¬ 
tic latitudes and dwarf galaxy surface brightnesses that 
are inaccessible with current detection methods. To do 
this, we simulated the RR Lyrae distributions of 13,200 
dwarfs, assuming numbers of RR Lyrae stars consistent 
with the RRab specific frequencies of observed dwarfs, 
and applied a simple FOF algorithm to select groups of 
two or more RR Lyrae with a 2d linking length of 500 

pc - . 

This approach can recover most dwarfs similar to those 
known that are more luminous than My = -3.2 mag (with 
a possible dip in detectability for My = -5 mag to -4 
mag, owing to observed specific frequency trends). Be¬ 
cause RR Lyrae can be distinguished from foreground 
stars near the Galactic plane, RR Lyrae could recover 
ultra-faint dwarfs similar to those known, but at low 
Galactic latitude. This approach also recovers nearly all 
dwarfs with surface brightnesses of 31 mag arcsec -2 and 
fainter that are more luminous than My = -6 mag, and 
will recover very low luminosity dwarfs (-4.5 ^ My <; 
-3.2) with high RR Lyrae specific frequency to similarly 
low surface brightnesses. Groups of RR Lyrae may thus 
provide the most effective means to te st the predicted ex- 
istence of stealth dBullock et alJl2010H and fossil galaxies 
(|Bovill fc Ricottill2011lh 

We used t he si mulated stellar halo catalogs of 
iLowing et al.1 (2015) to demonstrate that only modest 
contamination from groups of two or more field halo RR 
Lyrae is expected at d > 50 kpc, when excluding groups 
near major field halo structures. In particular we found 
(i) FOF searches with a 500 pc linking length are appro¬ 
priate to find dwarf galaxies with d > 100 kpc as groups 
of 2 more more RR Lyrae, and (ii) groups with as few 
as 2 RR Lyrae can also distinctly identify dwarf galaxies 
with 50 kpc < d < 100 kpc, however field contamina¬ 
tion should be limited at these nearer distances by more 
strictly excluding groups of only 2 RR Lyrae near larger- 
scale spatial clusterings of RR Lyrae. 

We searched several public RR Lyrae catalogs (CRTS 









































Dwarf Galaxy Searches with RRL 


11 


- CSS and MLS, LSQ, SEKBO, and Stripe 82) for rel¬ 
atively isolated groups of RR Lyrae to identify possible 
dwarf galaxy candidates with 50 kpc < d < 100 kpc. 
However, because these catalogs are very incomplete to 
RR Lyrae at d > 50 kpc, we can’t use this search to put 
quantitative limits on low surface brightness dwarfs. The 
Bootes I (66 kpc) and Sextans (86 kpc) dwarfs were found 
as groups of 4 and 2 RR Lyrae in the CRTS and LSQ 
catalogs, respectively, as expected given survey depths. 
This is far fewer than their known number of RR Lyrae 
because the CRTS and LSQ catalogs are severely in¬ 
complete at those distances. We associated a number 
of other relatively isolated groups of only 2 RR Lyrae 
with known halo structures, and identified two interest¬ 
ing groups that are not obviously associated with known 
structures. 

One of the most exciting regions of dwarf galaxy dis¬ 
covery space that RR Lyrae may open up is Galactic 
latitudes closer than ~25° to the plane, because the de¬ 
tectability of dwarfs with traditional resolve d star counts 
drops off rapidly at latitudes with b ^ 20° dWalsh et all 
120091 5. We caution that the simplified proof-of-concept 
calculations presented in this paper do not include the 
detection limiting effects of Galactic extinction, the in¬ 
creasing uncertainties in extinction (and therefore of RR 
Lyrae distance) with decreasing Galactic latitude, and 
the possible contamination of thick disk stars in RR 
Lyrae group finding at very low Galactic latitudes. Qual¬ 
itatively, Galactic extinction will decrease the effective 
RR Lyrae survey volume of each survey (especially at b 
20°) and distance uncertainties will soften the detectabil¬ 
ity of dwarfs. Even detection limits that are >2 magni¬ 
tudes shallower at the lowest Galactic latitudes, would 
represent enormous progress; More detailed simulations 
can explore these possibilities in the future. Including 
RR Lyrae periods in RR Lyrae group finding will help 
mitigate the impact of distance uncertainties and thick 
disk contamination on finding ultra-faint dwarfs at very 
low latitudes. Such dwarf s are populated by particularly 
long p eriod RR Lyrae ('e.g. lCatelanl2009l : lBoettcher et al.1 

MM. 

PanSTARRS 1 may be the most exciting existing 
dataset for discovering new MW dwarfs as groups of RR 
Lyrae stars. Its RR Lyrae survey volume could include 
~15% of the MW’s dwarf galaxies ('Section 12.21) . and it 
includes a large area at relatively low Galactic latitudes 
(unlike P/i/ZTF). If PanSTARRS 1 data produce an RR 
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Lyrae catalog that is reasonably complete to 100 kpc in 
areas of relatively low Galactic extinction, then groups 
of RR Lyrae at low latitude (or the lack thereof) could 
test the apparent significant spatial anisotropy in MW 
dwarf galaxies. The PS1 RR Lyrae catalog should also 
reveal candidates for extremely low surface brightness 
MW dwarf galaxies, if they do exist. 

Looking ahead, complimenting traditional search tech¬ 
niques with RR Lyrae in LSST may be the only way to 
obtain the most complete possible census of Milky Way 
dwarf galaxies - through the Galactic plane and to sur¬ 
face brightnesses of 31 mag arcsec^ 2 and fainter. LSST 
is the only planned survey that will be both wide-held 
and deep enough to grasp the majority of the MW’s halo 
with RR Lyrae stars. RR Lyrae identified in LSST’s 
wide-fast-deep survey will include ~45% of the MW’s 
predicted dwarf galaxy companions (even when account¬ 
ing for Galactic extinction). This quoted volume to be 
covered by the wide-fast-deep survey is a lower limit to 
LSST’s RR Lyrae survey volume, because ~10% of sur¬ 
vey time will be dedicated to thousands of additional 
square degrees of the Northern Ecliptic Spur, the Galac¬ 
tic plane, the Magellanic Clouds, and other special fields. 
RR Lyrae in the LSST footprint will be efficiently iden¬ 
tified to d ~ 600 kpc, enabling RR Lyrae-based maps of 
stellar structures and dwarf galaxies to well beyond the 
edge of the MW. 
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